Fluvial incision and tectonic uplift across the Himalayas of central Nepal

The pattern of fluvial incision across the Himalayas of central Nepal is estimated from the distribution of Holocene and Pleistocene terraces and from the geometry of modern channels along major rivers draining across the range. The terraces provide good constraints on incision rates across the Himalayan frontal folds (Sub-Himalaya or Siwaliks Hills) where rivers are forced to cut down into rising anticlines and have abandoned numerous strath terraces. Farther north and upstream, in the Lesser Himalaya, prominent fill terraces were deposited, probably during the late Pleistocene, and were subsequently incised. The amount of bedrock incision beneath the fill deposits is generally small, suggesting a slow rate of fluvial incision in the Lesser Himalaya. The terrace record is lost in the high range where the rivers are cutting steep gorges. To complement the terrace study, fluvial incision was also estimated from the modern channel geometries using an estimate of the shear stress exerted by the flowing water at the bottom of the channel as a proxy for river incision rate. This approach allows quantification of the effect of variations in channel slope, width, and discharge on the incision rate of a river; the determination of incision rates requires an additional lithological calibration. The two approaches are shown to yield consistent results when applied to the same reach or if incision profiles along nearby parallel reaches are compared. In the Sub-Himalaya, river incision is rapid, with values up to 10–15 mm/yr. It does not exceed a few millimeters per year in the Lesser Himalaya, and rises abruptly at the front of the high range to reach values of ∼4–8 mm/yr within a 50-km-wide zone that coincides with the position of the highest Himalayan peaks. Sediment yield derived from the measurement of suspended load in Himalayan rivers suggests that fluvial incision drives hillslope denudation of the landscape at the scale of the whole range. The observed pattern of erosion is found to closely mimic uplift as predicted by a mechanical model taking into account erosion and slip along the flat-ramp-flat geometry of the Main Himalayan Thrust fault. The morphology of the range reflects a dynamic equilibrium between present-day tectonics and surface processes. The sharp relief together with the high uplift rates in the Higher Himalaya reflects thrusting over the midcrustal ramp rather than the isostatic response to reincision of the Tibetan Plateau driven by late Cenozoic climate change, or late Miocene reactivation of the Main Central Thrust

Active folding of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal

We analyze geomorphic evidence of recent crustal deformation in the sub-Himalaya of central Nepal, south of the Kathmandu Basin. The Main Frontal Thrust fault (MFT), which marks the southern edge of the sub-Himalayan fold belt, is the only active structure in that area. Active fault bend folding at the MFT is quantified from structural geology and fluvial terraces along the Bagmati and Bakeya Rivers. Two major and two minor strath terraces are recognized and dated to be 9.2, 2.2, and 6.2, 3.7 calibrated (cal) kyr old, respectively. Rock uplift of up to 1.5 cm/yr is derived from river incision, accounting for sedimentation in the Gangetic plain and channel geometry changes. Rock uplift profiles are found to correlate with bedding dip angles, as expected in fault bend folding. It implies that thrusting along the MFT has absorbed 21 ± 1.5 mm/yr of N-S shortening on average over the Holocene period. The ±1.5 mm/yr defines the 68% confidence interval and accounts for uncertainties in age, elevation measurements, initial geometry of the deformed terraces, and seismic cycle. At the longitude of Kathmandu, localized thrusting along the Main Frontal Thrust fault must absorb most of the shortening across the Himalaya. By contrast, microseismicity and geodetic monitoring over the last decade suggest that interseismic strain is accumulating beneath the High Himalaya, 50–100 km north of the active fold zone, where the Main Himalayan Thrust (MHT) fault roots into a ductile décollement beneath southern Tibet. In the interseismic period the MHT is locked, and elastic deformation accumulates until being released by large (M_w > 8) earthquakes. These earthquakes break the MHT up to the near surface at the front of the Himalayan foothills and result in incremental activation of the MFT

Stress transfer and strain rate variations during the seismic cycle

The balance of forces implies stress transfers during the seismic cycle between the elastobrittle upper crust and the viscoelastic lower crust. This could induce observable time variations of crustal straining in the interseismic period. We simulate these variations using a one-dimensional system of springs, sliders, and dashpot loaded by a constant force. The seismogenic zone and the zone of afterslip below are modeled from rate-and-state friction. The ductile deeper fault zone is modeled from a viscous slider with Newtonian viscosity ν. The force per unit length, F, must exceed a critical value F_c to overcome friction resistance of the fault system. This simple system produces periodic earthquakes. The recurrence period, T_(cycle), and the duration of the postseismic relaxation phase, which is driven dominantly by afterslip, then both scale linearly with ν. Between two earthquakes, interseismic strain buildup across the whole system is nonstationary with the convergence rates V_i, just after each earthquake, being systematically higher than the value V_f at the end of the interseismic period. We show that V_i/V_f is an exponential function of α = T_(cycle)/T_M ∝ Δτ/(F – F_c ) ∝ Δτ/(νV_ 0), where Δτ is the coseismic stress drop and V_0 is the long-term fault slip rate. It follows that departure from stationary strain buildup is higher if the contribution of viscous forces to the force balance is small compared to the coseismic stress drop (due to a low viscosity or low convergence rate, for example). This simple model is meant to show that the far-field deformation rate in the interseismic period, which can be determined from geodetic measurements, might not necessarily be uniform and equal to the long-term geologic rate

Modeling mountain building and the seismic cycle in the Himalaya of Nepal

A host of information is now available regarding the geological and thermal structure as well as deformation rate across the Himalaya of central Nepal. These data are reconciled in a two-dimensional mechanical model that incorporates the rheological layering of the crust which depends on the local temperature and surface processes. Over geological timescale (5 Ma) the ∼20 mm/yr estimated shortening rate across the range is accommodated by localized thrust faulting along the Main Himalayan Thrust fault (MHT). The MHT reaches the surface along the foothills, where it is called the Main Frontal Thrust fault (MFT). The MHT flattens beneath the Lesser Himalaya and forms a midcrustal ramp at the front of the Higher Himalaya, consistent with the river incision and the anticlinal structure of the Lesser Himalaya. Farther northward the MHT roots into a subhorizontal shear zone that coincides with a midcrustal seismic reflector. Aseismic slip along this shear zone is accommodated in the interseismic period by elastic straining of the upper crust, increasing the Coulomb stress beneath the front of the Higher Himalaya, where most of the microseismic activity clusters. Negligible deformation of the hanging wall requires a low apparent friction coefficient (μ) less than ∼0.3 on the flat portion of the MHT. On the ramp, μ might be as high as 0.6. Sensitivity tests show that a rather compliant, quartz-rich rheology and a high radioactive heat production in the upper crust of ∼2.5 μW/m^3 is required. Erosion affects the thermal structure and interplays with crustal deformation. A dynamic equilibrium is obtained in which erosion balances tectonic uplift maintaining steady state thermal structure, topography, and deformation field. Using a linear diffusion model of erosion, we constrain the value of the mass diffusivity coefficient to 0.5–1.6×l0^4 m^2/yr. This study demonstrates that the data are internally consistent and compatible with current understanding of the mechanics of crustal deformation and highlight the role of viscous flow in the lower crust and of surface erosion in orogeny processes on the long term as well as during interseismic period

On the growth of normal faults and the existence of flats and ramps along the El Asnam active fold and thrust system

The combination of detailed topographic leveling on the southwest segment of the El Asnam thrust fault with existing seismic and geologic data implies that the geometry of this fault involves shallow dipping flats and steep ramps. The fault appears to be growing along strike toward the southwest end, where the main shock initiated in 1980. From a depth of about 10 km, the main thrust appears to ramp to the basement-Cenozoic cover interface on a plane striking N40°E and dipping 50°–55° to the northwest. Along the southwest segment where folding has not yet developed, the thrust continues steeply through the Cenozoic cover to the near surface where it flattens, causing normal faulting. Along the central and northeast segments, which display a more evolved fold structure, the deep thrust probably flattens at a depth of 5–6 km, into a decollement along the Cenozoic-Jurassic interface before ramping to the surface. The Sara El Marouf and Kef El Mes anticlines have thus formed as fault propagation folds. Normal faults at Beni Rached probably branch with the thrust to maintain kinematic compatibility between the deep ramp and decollement. The greater separation (∼7 km) between the normal faults at Beni Rached and the thrust where it crosses Oued Cheliff than along the southwest segment (∼1 km) reflects the greater depth of the ramp to flat bend. We infer that the September 9, 1954, earthquake activated only the central deep segment of the main thrust together with the Beni Rached normal faults, while that of October 10, 1980, activated the whole system of flat decollements, ramp thrusts and compatibility normal faults. Further complexities of the faulting in map view are related to changes of strike of the thrust (in particular north of Oued Cheliff)

Erosion as a driving mechanism of intracontinental mountain growth

In nature, mountains can grow and remain as localized tectonic features over long periods of time (> 10 m.y.). By contrast, according to current knowledge of lithospheric rheology and neglecting surface processes, any intracontinental range with a width that exceeds that which can be supported by the strength of the lithosphere should collapse within a few tens of millions of years. For example, assuming a quartz-dominated crustal rheology, the relief of a range initially 3 km high and 300–400 km wide is reduced by half in about 15 m.y. as a result of lateral spreading of its crustal root. We suggest that surface processes might actually prevent such a “subsurface collapse.” Removal of material from topographic heights and deposition in the foreland oppose spreading of the crustal root and could eventually drive a net influx of material toward the orogeny. We performed a set of numerical experiments in order to validate this hypothesis. A section of a lithosphere, with a brittle-elasto-ductile rheology, initially loaded by a mountain range is submitted to horizontal shortening and to surface processes. If erosion is intense, material is removed more rapidly than it can be supplied by crustal thickening below the range, and the topography is rapidly smoothed. For example, a feature 3 km high and 300–400 km wide is halved in height in about 15 m.y. for an erosion coefficient k = 10^3 m^2/yr (the erosion rate is of the order of a few 0.1 mm/yr). This regime might be called “erosional collapse.” If erosion is not active enough, the crustal root spreads out laterally and “subsurface collapse” occurs. In the third intermediate regime, removal of the material by erosion is dynamically compensated by isostatic rebound and inward flow in the lower crust so that the range can grow. In this “mountain growth” regime the range evolves toward a characteristic graded shape that primarily depends on the erosion law. The erosion rate may be high (e.g., 0.5–0.9 mm/yr), close to the rate of tectonic uplift (e.g., 0.7–1.1 mm/yr), and few times higher than the rate of topographic uplift (0.15–0.2 mm/yr). These experiments show that surface processes can favor localized crustal shortening and participate in the development of an intracontinental mountain. Surface processes must therefore be taken into account in the interpretation and modeling of long-term deformation of continental lithosphere. Conversely, the mechanical response of the lithosphere must be accounted for when large-scale topographic features are interpreted and modeled in terms of geomorphologic processes

Geodetic displacements and aftershocks following the 2001 M_w = 8.4 Peru earthquake: Implications for the mechanics of the earthquake cycle along subduction zones

We analyzed aftershocks and postseismic deformation recorded by the continuous GPS station AREQ following the M_w = 8.4, 23 June 2001 Peru earthquake. This station moved by 50 cm trenchward, in a N235°E direction during the coseismic phase, and continued to move in the same direction for an additional 15 cm over the next 2 years. We compare observations with the prediction of a simple one-dimensional (1-D) system of springs, sliders, and dashpot loaded by a constant force, meant to simulate stress transfer during the seismic cycle. The model incorporates a seismogenic fault zone, obeying rate-weakening friction, a zone of deep afterslip, the brittle creep fault zone (BCFZ) obeying rate-strengthening friction, and a zone of viscous flow at depth, the ductile fault zone (DFZ). This simple model captures the main features of the temporal evolution of seismicity and deformation. Our results imply that crustal strain associated with stress accumulation during the interseismic period is probably not stationary over most of the interseismic period. The BCFZ appears to control the early postseismic response (afterslip and aftershocks), although an immediate increase, by a factor of about 1.77, of ductile shear rate is required, placing constraints on the effective viscosity of the DFZ. Following a large subduction earthquake, displacement of inland sites is trenchward in the early phase of the seismic cycle and reverse to landward after a time t i for which an analytical expression is given. This study adds support to the view that the decay rate of aftershocks may be controlled by reloading due to deep afterslip. Given the ratio of preseismic to postseismic viscous creep, we deduce that frictional stresses along the subduction interface account for probably 70% of the force transmitted along the plate interface

Tectonic evolution of a continental collision zone: A thermomechanical numerical model

We model evolution of a continent-continent collision and draw some parallels with the tectonic evolution of the Himalaya. We use a large-scale visco-plasto-elastic thermomechanical model that has a free upper surface, accounts for erosion and deposition and allows for all modes of lithospheric deformation. For quartz/olivine rheology and 60 mm/yr convergence rate, the continental subduction is stable, and the model predicts three distinct phases. During the phase 1 (120 km or 6% of shortening), deformation is characterized by back thrusting around the suture zone. Some amount of delaminated lower crust accumulates at depth. During phase 2 (120 km–420 km or 6%–22% of shortening), this crustal root is exhumed (medium- to high-grade rocks) along a newly formed major thrust fault. This stage bears similarities with the period of coeval activity of the Main Central thrust and of the South Tibetan Detachment between 20–16 Myr ago. During phase 3 (>420 km or 22% of shortening), the crust is scraped off from the mantle lithosphere and is incorporated into large crustal wedge. Deformation is localized around frontal thrust faults. This kinematics should produce only low- to medium-grade exhumation. This stage might be compared with the tectonics that has prevailed in the Himalaya over the last 15 Myr allowing for the formation of the Lesser Himalaya. The experiment is conducted at constant convergence rate, which implies increasing compressive force. Considering that this force is constant in nature, this result may be equivalent to a slowing down of the convergence rate as was observed during the India-Asia collision

Craters as sand traps: Dynamics, history, and morphology of modern sand transport in an active Martian dune field

Aeolian transport of sand is abundant on modern-day Mars, as revealed by remote sensing measurements of the motion of dunes, and of the meter-scale ripples that mantle them. We study a large-scale natural sand trap within the Meroe Patera dune field: a 1.8-km diameter crater which features a dune-free “shadow” in its lee. We compare the volume of sand trapped within this crater to the sand volume that would be expected to cover the area of the crater and its dune-free shadow behind it if the crater were not present. We find that the crater holds less sand than this “missing” volume would predict, implying that sand escapes from the crater over time. Modern day imagery shows an apparent lack of sand escaping from the Meroe crater, however, suggesting that changes in the wind regime at the site may have allowed sand to escape in the past. The persistence of an altered dune morphology all the way to the far downwind edge of the dune field suggests consistent wind conditions over the time of the crater-dune field interaction

Investigation of the relationships between basin morphology, tectonic uplift, and denudation from the study of an active fold belt in the Siwalik Hills, central Nepal

The present study investigates correlations between an extensive range of geomorphic properties that can be estimated from a digital elevation model and the uplift rate on geological timescales. The analysis focuses on an area in the Siwalik Hills (central Nepal), where lithology and climate can be considered as uniform. This area undergoes rapid tectonic uplift at rates of up to 15 mm yr^(−1), which are derived from the geometric pattern of a fault-bend model of fold growth. The selected geomorphic properties can be divided in two categories, depending on whether or not the vertical dimension is taken into account. None of the planar properties are significantly correlated to uplift rate, unlike those that include the vertical dimension, such as the mean elevation of basins, hypsometric curve, and hypsometric integral, and relief defined by the amplitude factor of length scaling analysis. Correlation between relief and uplift rate is observed for all length scales of topography shorter than 600 m, which suggests that all orders of the streams are able to adjust to the tectonic signal. Simple mass balance considerations imply that the average elevation is only 10% of surface uplift, suggesting that a dynamic equilibrium has been reached quite rapidly. Using a simple two-process model for erosion, we find that fairly high diffusion coefficients (order of 10 m^2 yr^(−1)) and efficient transport of the material by rivers are required. This unusually high value for mass diffusivity at small length scales may be obtained by either a very efficient linear diffusion or by landsliding. Actually, both processes may be active, which appears likely given the nature of the unconsolidated substratum and the favorable climatic conditions. Local relief in the study area may therefore be used to predict either uplift or denudation, but the prediction is calibrated only for that specific climatic and lithologic conditions and cannot be systematically applied to other contexts